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(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a parallel matrix 
processing method suitable for a shared memory type 
scalar computer. 

SOLUTION: A matrix to be LU decomposed is divided 
into a block D at a diagonal part and a row-direction 
block below the D, for example, L1 to L3. D+L1, D+L2 
and D +L3 are allocated to three processors, 
respectively, to perform operation in parallel. The block 
U is updated by the method of LU decomposition, and 
further C1 to C3 are updated using L1 to L3 and U. The 
processing is conducted for the inside of the block 
sequentially decreased in size, thereby lastly leaving the 
part corresponding to the D. The D is LU decomposed 
by one processor to end LU decomposition for the whole 
matrix. 




LEGAL STATUS 

[Date of request for examination] 19.12.2003 

[Date of sending the examiner s decision of 
rejection] 

[Kind of final disposal of application other than 
the examiner s decision of rejection or 
application converted registration] 

[Date of final disposal for application] 

[Patent number] 3639206 

[Date of registration] 21 .01 .2005 

[Number of appeal against examiner's decision 

of rejection] 

[Date of requesting appeal against examiner s 
decision of rejection] 

[Date of extinction of right] 



http://www19.ipdl.inpit.go.jp/PA1/result/detail/main/wAAAVLaypIDA414163246P... 2007/03/31 



(ig>B*n*ff* (jp) 02) & ^ ^ ^ (a) <n>»afaM&H## 

#^2002-163246 
(P 2002— 163246 A) 

(43) ftW H ¥J3U4*E 6 7 B (2002. 6. 7) 

(5i)intci. 7 wwe* FI f-?a-r(##) 

G 0 6 F 17/12 G 0 6 F 17/12 5 B 0 4 5 

15/16 6 1 0 15/16 610Z 5B056 



SSEM* »^3S©»7 OL (£12J€) 



(21) fflH#^ 

(22) ffiljB 



)-358232( P2000-358232) 



JF«12¥11H24B (2000. 11. 24) 



(71) fflKA 000005223 

i# 

(72) 5§W# *B M 

1# S±a#:££*fcl*9 
(74)ftSA 100074099 

#s± mt (m«) 

F^-A(#55) 5B045 AA07 DD02 DD12 GG11 
5B056 AA04 BB02 



64) mmco&mi p&*&mx%?&mt&m\z&vz>mmMsmJiiii. nzftmim 



(57) [Sift] 

Di:Dfl>TfiiII=fcii5iJ^lSjro^P'V^. CQAIf. L1~ 

iz. d+li, d+l2. D + L3$iiiyay. 

jSlCfco-CMfcrU L1~L3iUtSffl^tC 
St -5 c <t ir «fc y ^5iJ±«:ro LU»SS8x§. 




L2 



U 



Ul ! U2 | U3 



C2 



(2) 



2002—1 63246 



[tt#Sil] **0^P4sv«-*S?a-^*1»O**jC 
£'hfr5>J ? § ?P ? ftXT- v -Jt . 

Kaiaro^p-fe^-y-^vi-^afeS'iiZx **i-eMr-r 

[»#H2] gtfj$iJ5SIU±. ^5iJ©LU»St?fe*Ct 
£4#St<fc-*-§W#JSl lzE«<7)Eii&f*,, 

[»#js3] ffjgasaio^p-b^-y-^i? 3. -juA<wr* 

/J^5"J^P vVO)*?— $ frb-eft-eh-tftftfy F<D«*§ 

KbfTtw hansffl*^Ka?aro^ , p-b-y* a Evn.-;n- 

WfcTJUa* 'J XAI- J: oTIH&fr^JiDMSfrSmv WE 
■9- ; Evi-;ut?LU»<g-r-&^i:|-<fcy. ttfr*J±i*i:: 

•Di*-cLu»»**7-r*ci:*iMit-r*i»*a2i= 

EttttEIMKtt. 

[St#JS5] W&WSimt. fi$iJ(DP UX*-#fl?fci> 
1 1cKKOES&t*. 

[Bl*«6] BtrEPl/X4r-##fcSlM*£»=iUX* 
-5MII*. ffJEfx5<Ja>*t-ffl]fr &!WBWfc7JU =f 'J XAI= 
J:oT«J*fTW«HE**fTlV flfEfr^JOrt. SttKX 

»»-*-4Ctl=«fcy. »fr9J£tt;l=oi\-ci-U#»£5S 
T"T * - <t VMHt-t -511*51 2 lzE«<BE»»H*. 

fr*"J ? p v * a =flfraW# * frEa»ro zfa -tr -9-=E 
i?i-;KO»<D2fgOD»05^p^^|r»fiJL. 
*ifcH^fr5i]S?^0)^P yf$20fOfi^fC, 
#^P -b v 5? a.— P — * * ^ 'J ttlcfetfl U 
?a-ky1*-«9a.— JM=fy*rtt« = t Zft&t-f 
SE*#i3 5 izEttaBftltlt. 
[fflOMUMU] 



[000 1] 
[0 0 0 2] 

[«£*<*>&»}] it£l*a«a*tl-JHlfcJ:o-C«<» 

[0003] lii«*§att. «**«f 

JHcSf U<Jjt*iL^3J:5l=Ea-r-6Ci:)!i<a**. C 
ZV. LU»S (LU Factorization) il^3^acl=<fc*i 
I*. «S!t$S-rfr5"J*±H#lff5iJ (upper-triangular ra 
atrix ) bT=.fk?i&\ (lower-triangular matrix ) iz 
ttM-fZzkizk-DX. iiftljfcaas*©****!)*. 

t&zt&msL 1 **a*«>****fc*>i=Mefc»a 

UX*-»S (Cholesky Factorization) <t^5ff5»J<» 
a ) SfrMail£ 1 «»«0|HI 

n«5ija>5iitiafe^asc©»a?i=HLr, h*jfcM 
•+jwc©iijci**a*i*. ^pv^ftLtntao 

^*^JMaf=^P«y^*LU5MB (1) LfcSL fcffcL 
fcfr^f h)U$mtefc^P-v^$SESrUT (2) 

jE*^frw*H« o) -r«aa«aua-r. 

[0004] (1) ©5B31(4, — o©^P-bv9- 

^fSltl 2 («**^l*fj : «8MTM«>fMSfc*lM*5<J 

«**-r) a*<DJta«/h*^ai:L-cufc. cwiss 

(2) (3) <OfiP^(DS«rt,i(ei 2aflE©ffW3t* 
[0 00 5] fttaXKO**^ (3) ©l+JT-C, 

^'jiaxa^a^Ji+JMa (smp) -ei*. 
[0006] -r^*.. o) ©i+*i*. ftjua-e* 

tlTLN*) JP-KLT. ISSfSSSX HT-r*i^3 

a*i=jt^?**<ttyttfii«5i*uitttti.*. 

[0007] zotzih. -Jay91&£±£<-rz>&&:t>< 
?P'y?*S£*£ ^©^P-^(DL 

b) jEffi*t^fi-9l©aSt 1 *7Sa3t<DS£ 

iEflS*j«;ft5ij©aa i *7jasro»^icgsi lt. tha 
5Mk» % 'J anjut+ffa-eii^ftM^p v f t o y c i 
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/h*<-e#. at»<fc«i**Ka<>4ct3Ei<Rr»-cfeofc. 

[000 8] 

Slxa^aMMMMTtM** LU»Sfc^>i>iizj ux^r- 
xh*y** *#**yic7$-feX-*-ftfcJ6a>:aX htf 

[0 0 0 9] ##gB.S<Dg®l*. **M^U8!X*^tU* 
[0 0 10] 

»lt*^n^^ftXT-^^t. tt/hfrM^Ps/£<0l*k 

[0 0 1 1] **©^P-b^-»*»a 
— ;u*<**ff*t*-r*p— *;u>^y««£*rSbl^PIJB 

[0 0 12] 

[ooi 3] **>*yax*5ifc«tt««i*. 

?P-bV»1 0-1. 10-2. ■ • ' 10-ntf2S 
^*r>a^'J 13-K 13-2. ■ ■ ■ 13-n 
t^LtUMMi 2(cma*n^>o **P-fe*1M 
0 - 1 . 10-2. ■ • • 10-nlt ^-(DrtSfffo&l^ 
I*. 2****;^i>»*"J13-K 13-2. • ■ ■ 
i3-n*«?Ptv»«i:i r>a>^'J *<K 
It6*i4. «fc» #?p^*10-K 10-2, • 
■ • 1 0-nlC#:&£fcoTt^*^y^&S?i— 1 
-1.11-2. ■ • • 1 1 - n l±. ffl5$S£-J|?3 1 2 £ 
^Lr^P-b^-9-1 0-1. 10-2. ■ • • 10-n 
^T^-trXplflt^-Dt^o 10 

-2. • ■ ■ 1 0-n#T-*»a*ff3»*l=tt* * 

t\ /^'J^yjL-jH 1-1.11-2. • ■ -11 



-n«*& 1 O0)^P'fe^-9'A<a^-r*T r -^*2^-V 
«r>a>^'J13-1, 13-2. • • -13-n \Z& 

$JL. Sic. a^^^^^'J^&ftWfitW 

[0 0 14] Mi^5fc, iX^^i/^^'J^ 
&2**^!y*>i^ i EyiC«ka-T-*^»tt4tl. 2* 
+ t7va^^U P^0>^— *tH±X$kmi>%k*>Z>£. * 
t'J^va-JM 1-1.11-2. • • ■ 1 1-n0 

«■ :?p -fe ^ l fc^o) x—* £ > ^ y * v a.— ;u 
^bf-^sa^tig, ^^PHsv-s-^a^o^-f s 

LfctfoT. Z<7)i:#lCfi. ^TO^P-b^-y-^y^E'J 

^+*-<D«i31(^l^^£<t^Ci:£/<y Tl^^ (Barrier S 
ynchronization) £5&£><tL^do 

[0015] S2st;0 3f*. ^szwommmmizft'D 
i*. laa-r^sfrWflowiCB-efcy. B3i*. aisma 

*nj6^®iwfci>rii. ^p-tr^-tr^sas^aa-r*^ 
#£ftmt-?z>o Triztt. 0 2!zfci*r. «*-ci** 

L 1 ~~L3<&gfl#l*. ?P'^B^h*<Lt. no 
^P-tr^-9- (PE. &£IM*. XU*yK) -CfliSLri> 

■tZt&\Z^ L 1 -L3£-ett-m»JC0XLx-; KlcSiy 

[0016] ^W^ : EyMx*v^Jft»:^«.^P4? 

y it IC 1 2 2fc V *y V J. ifiVt t> o T I > £ o 

5 d tfxttitt £ 1 1 ^ m-r Ji -e fc £ o 
[0017] T : -*fi(7)**^p ( aa£^»PE-e«<ii 
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46 £ JtSf -5 z: t tf&m t * . 

[0 0 18] COfctolC 0 3IC^**Vi>*3l-. 

p-b yH-na— a;u=i*****«LT L 2 v ->i 
&5"]|c£^£,h#i&gfc:?P-;/?*ta8l5#DI;Ju 

[0 0 19] b:^P-y^«. k : «-^P-b-y-9-A^fi 
"T -55IJ-?P v +> <0 1 S <D*£ £ C?P ? tStf/jNfr 

5"J^p-v^<D5iJ^isi-efc*©-e. '^'n&l-?a-y<?o?io) 

ft. -f£*>*>. 02l=j3ltSL 1 ~L3OT^P-y5'(D^ 
tf-<0fj<08a) i:Lfci:#. bx (b + k) xs~8Mb 

[0020] fit, tt3S« (IJH^yi/a/t'J) 

d5*-f£* (fr5"JA<*#<) , **il=tt 

^T^P-b-v9-aA<^*^fc«>. 365"J«iSlRlltaJ5iJ^P 
^^<0^P^^tSA^J>*<*-5i:#li, ?D?>IS» 
9IL,-C^gJS:^P^-j7*«i«LT3!)^. ff^Pv^fflS 
•r & fTaM«>X«r£ttHI=fT 5 . 

[0021] mi-, ft^p-t-y-y- (*xu»yK) -eroft 

u«-S£tr5^*£*Jffl-r-S>w£-e. *+-r>a±fflf 
— £ <*Jffl-f -5. 

[0 0 2 2] H4li, *lt&fl2JgOTLUtt8?cDfliJia)3iE 
H^Tjk-rHSt^va— =f-*t— t-X'&Zo *-*\ X-rv^S 

ff5u©*#*) ^e>:?pv$*B£*s-f x^ 

£L. K-C«ia-r*^P-yJ7 (Sl2(7JDaT/L 

i) *f***icpe— r-5o fir. xf->?s3i:fc 

nr. #xu>v K-ctf#? f*£*5SL. -€-ro+-efl)ftX 

O 

[0 0 2 3] X^r^y^S 4lCfclNTIi. flHfcMfcibyrt* 

^Ui (H2»W) SLU»»0>T^UXAlCtt-3 
TJE*r*"So fit, 6lz£LvC» 

K-eat«ic. ?p^c i (B2#n) t?p^n 



tZJUytU (U i *«^*:b1±fc*a> : H2#i) (0 
«r*M»r-r*o -fLt, X-r^S2lcRy % &CD:?P 
(Ci^bftS?P^i?) ^8(DSatLU» 
HW*. fit, *KlC/Jv*<3S4*fflS^P^^^« 

aiz^p^^Dic»jt-r*»»iD^i=ay. noxu 

[0 0 2 4] H5-H1 or*. **lfi^J8<DLU»«a) 

^atiUBHtHMBt-SHtii. ccoi*. 204 

8 x 2 0 4 8 0)f7$lJ£4Xl^ K"CLU»IBt4»«« 
[0 0 2 5] Sf, 2048x2048(DWJS^n7 

^fi2 5 6tfLU^«t^t(?)i:no fr?n?94>B 

»iti±. Hsrc^-rtfey-e**. fir, 4oocp 
u (xi/*; K) K-ea 

8fcfc^Ii& ( (256 + 448) X256: 8MB (L2 
*+visa.to**m *y/h*l^««) *«M*U H6 
fcaM-J^lC #XU? K(7)#ffl«lCD 1 +L 1 „ D1 
+ L2, D1+L3, D1+L4?Je-«, 

[0 0 2 6] ^P^^«(*«xtf. J2lT<D<fc3l:: 

LTSfe^o Fp^S0>*£* (ffWOD*** : frM0)A 
ft) £n. K<Dtt*#TH READOUT. 

[0 0 2 7] 

[»1] 

l/Cn* /# THREAD) x — - nb 
i 100 

[0 0 2 8] ttSl^T. 
nb^51 2ftb, :?P -/ ^7 if§= 5 1 2 
nb^2 5 64*b, ?P?£M=2 5 6 
nb^1284b, ^Dy*«=128 
fftJUJk ^P'^B=64 

[0 0 2 9] tft^, LU^S03Xhlt 2n3/ 

T, 4ftfl)3XhS#THREADT*fML, 

IX^KW^p*^ ffiSS t * J: 5 

[0 0 3 0] zzv. mfR£tt}\i*tzto* ttMflsLftt^ 
»$0)LU»»<7)T;U=rUXA*H7|C^-ro S7lCfc 
L^ri*. LT = D 1 + L 1 + L2 + L3 + L4 (D§fl#£ 
LU^St^T^UXASfvLtt^p LTlt 2 0 
4 8X 2 5 6©?Py*ift'3tl^. 

[0031] (i ) a)»»rcfe^-c. h* 
*s-r^>o ibiks iiLTO)fit*y, 256 

2048-Cifc^ j J 1=1*. e#^Hfl)«t«fiM 
[0 0 3 2] fit, (2) (DfflJ^lCfcl^r. ASAS 
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j j i SCDfT^^— * £ j j 

frtf)^-* tAtiWx-Z>o 'X\z^ (3) (DSP^-e. W i 

[0033] ^0 (1) - (3) £ i #1 ~ iblks izfr 

fcoT«ysL3MW*. cct?. LT(D§*-r?feS, 

leng*t*oi*:*<tt«i:. ZtitommitL 2*-V v 

H8lc^-Ti:j3y-e** 0 
[0 0 3 4] fcfc* H8lzj3l*TI*. LTilt * 
pivot (4). GPIVOT. ROW (iblks ) \t. &m®\z 
»ttS*i*5*—*-T?**. (4) lcfct>-c. #x 

\sv K-Ctf#v b£3*£ 0 fit, (5) trS^fil^iS 
$i]pivot(4)(DXl^ Ki§a)SSI^fitttl>o (5) <0 
fftlc. /<'J7HH*Sot, (6) fcfcl^T. S^tf/K 
* h*»oxu* K»*£GPIVOTfc»IA-r*. *LT\ 

(6) <&&t*. wtf/^UTHJHsa*. (7) rc 

(8) icfe^tl*. GPIVOT*<OC0<!:#I*. **b?tf^ K 
l*D10>rtl=fc*fr* A*t»Jt^FK"C*y* P— 
IC A*ltf X 4 » GP I VOT## T H R E A D \Z*$ L I* k & . 
-r*t)*». GPIV0TA<O*y**^»dlCtt. t*tf#f 

HU£A*lffx.£o fit, (9) % (1 O) iCfet^ 
T. LU^»<Dfc^0);1ig:^fT5o fit, ±SB (4) 

- do) ^mmrt^nv^o±xa>Wz-Di^x^7 

5o *rft*3%. 1 ~ ibiks *-ccd i iz-3L\-cdnasi»y 

[0 0 3 5] :Ct\ LTKOLU^S^'JT 
l^»!£l*&o K©D 1 <Dffl5#(£L L 

£uulcLU##S*iTl**. fit, ftXU^K-eu 
K-LL-1U1, U2<-LL"1U2. U3<-LL^U 
3. U4<-LL-1U4^U^ K-CatMfZtW-*"*. 
C(Dt+S:(D^r% DK LK L2, L3, L4£P- 
*;b«^61¥«Alcae— /<v^L, /<'J7WIJH*tt 
£o Hf^ C1*-C1-L1XU, C2<-C2-L2x 
U. C3<-C3-L3xUs C4<-C4-L4xUSt 

^p^/^«iblk8*«/hLfc»»SHi:«fc^l=»«IL. 

g!9lC^-T<t D1, LK L2. L3. L4<75^ 

I^5c C(DJ:5lc. aaeitySLTl*< 



Xl/«> K"CLU#»LT«3S*»T-r*. 
[0 0 3 7] W« ±IBffl3lT7(±. LTijLU^jgt 

z>tztb&t%&}teLuitf»£mmLx^z>o *fc. c i<- 

[0 0 3 8] HI 014. S«WLU»»7^U^ 

^-f>LUt Lt4^b^o LU<DIfc££Sl£. 

D 1 + L i &«M) , k (LTi 0)l#.5tS<D*L$ 
*) . iblks C?P^«S) (DttlC. LU#fi££*&tf>3 
ttS^^-fist. LU^S^^T5raT?fc*nwidr*fe^c 
[0 0 3 9] -th^U— T>0)fitO-e. nwidrf><8. 

■tat)*). LU5MB*fr3«3B<8-efc*3&i5*tflBf-r 
*« YES0)4£lwti:* LT i (ist : k. ist : ist 
+nwid-1) *3M«lcLU$MiM-4. z^X. ist : k 
kmoBORit. Sadist rf^kST-cDLT i 
mx % ist : ist +nwid- 1 kl^otol*. SE»J&<istfr 
fcist +nwid- 1 tt?(DLT i &7ji'?M&X&Z>o J*IT 

[0 0 4 0] LTi (DLUflJBIZfclvCI*. ±IB (4) 

-do) (Dmm&'noo fceu *7<DA*Lfixatf* 

I*. g£ iblks LTi ( i . 1 : iblks ) fA^lf 
7L%> 0 ±EttW*<NOa>»«l=li\ LU^ft 

5 ignwi d ^ 2 T*fO o fefllCO L U CD L U ^»CD-9- ^ 

;u-^>^w«MicP?t;aiLT^f5o trs^ 

^>*i»tfart". Z(D;u— LTi (is 
t : ist +nwid/2-1. ist +nwid/2 : ist +nwi 
d) *S«r-r*o Hl^. LT i (ist : ist +nwid/2 
-List : ist +nwid/2- 1 ) (DTHftfi^JL L $■ 
SfflLt, LL-1jCi lc£fr£*M*THSrf £o ^ 

[*. 

LT i (ist +nwid/2 : k % ist +nwid/2 : ist + 
nwid) =LT i (ist +nwid/2 : k. ist +nwid/ 
2 : ist +nwid) -LTi (ist H-nwid/2 : k. ist 
: ist +nwid/2- 1 ) xLTi (ist : ist +nwid 
/2 — 1. ist +nwid/2 : ist +nwid) 

*LT. /<U7Wlt»^ LU 

b) jEtt»»fT«a>iisc i **asco)»a 

■ 1 1 —S1 1 3 1*. E«»Wff>Ja)»*l=aUX*— » 
[004 1] ^fi5ij^a^:i ^Sa«tB»l^. fr*U£ 

D* LK L2. L3l^9lLT. *XU*y KlC*fft^ 
SlDfcMrtW^P^»»L1, L2. L3SM 
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P'^s^j3^+-^«t^o ate. z0>«^. tr 

fiJJBLT. /hTHft^J (C1-C63d>b^) £E« 

£L*z<t#. 2x#Ticmi:?n*2«ic#it* era 

;fc>*>. C1tC6, C2<hC5. C 3 <h C 4 

2x#T + l- iilO^P 7**Hft4^itft« 

[0 0 4 2] «S*iL-C^4fiW*<jE«»»fTWe*« 
H2 0UlZ»Ifr*»»l*. L1+L2+L3* 
fettiM^n^tDlEltLT fctt«D-Cl*5<D"C CCDS 
*|Z(±. LT 3S**ff5^S^*<, LI. L 

2. L3*-t*L-P*tK«LT3bf--r*lli«feU^ 

[0 0 4 3] H1 1 31*. S«»=l UX^r-^S^SaSCD 
«fra>tt31*IttB**B-C**o (1) IwJSl* 

T. frM<D-#3£a>^a^^*a L/X*-#»L. 30 
0^*3 1 iC^bf— "f£o fit, 3 1 a>Ta>J&Kn? 

(2) icte^T. M-r«M«EMI«2«lcLr3 2Q« 
#£3 3lz=jfcf— L. 3 3<DT<0jSJ&T-H£;ft£:gB#£ 
*M»»#*JBlvC\ E$rlT*o Its (3) 
T*. 3 4 0>aJ#£3 5ICZJ t°— L. 3 5 (DT0)jjklt^ 

(4) ICfeLNT. 3 6-e^^tL^>fiP^*3 7IC=je— 
L. 3 7<DT0)ffl»*«tt»»*ffl^T3S«T-r*o E 

(5) iCtel^T. 3 8 £ 3 9 lea If— L. 3 9 <DT 
a)Altl?w**iTl^*»»*l»litt»*ffl^TE«rL. 

(6) leftist. 4 4 1 iCZItf— L. MMFfteiH 

ivc. 4ia>T0>»#*E«*-*. EI-. (7) icteU 

T. 4 2 0>»#*4 3|C=ltf-L. fMi»#*ffll*T4 
UX*-»»*S»Mfc«yfiLafflLa*<e>. Attn 

[0 0 4 4] B1 4-H1 6 I*. ir^l^X^-^SO) 
I*. RW*«BMb**fctf>4XU* Kt*oT»31tfT 

[0 0 4 5] *-T. #Xb7 KI=H1 4fC^£;K£Dx 

&tfL i ««tnaiiMLT i icne-t^o fit. 
lt^ldlt ##-r LDJ »»i*n»Ma^ 

^T*fr5o fit, D L i ^J2lT0)ttjr-C«tatWI=3 
[0046] DLi<-DXLiK CZT?. Dl*Dx<0 

lc. t,i:0!)fil«lc=je— fir, ::-c, 



8. C2<tC7. C3tC6. C4iC5$^7i:L 

• xb7 K i -ei*. 

C 1 «-C 1 -L 1 1 x D L 1 1 
C8<-C8 — L4 2 x D L 4 2 

• K2-eii. 

C 2<-C 2-L1 2XDL1 2 
C7«-C7-L4 1 XDL4 1 
■ Xb'V K3"Cli. 
C 3«-C 3 — L 2 1 x D L 2 1 
C 6<-C6-L3 2 x D L 3 2 

C4<-C4-L22XDL22 
C 5<-C 5 -L 3 1 x D L 3 1 

4:1*5***159. 

[0047] -tux. zc*T?a>**#ilfrofcB*jiii-e 
/< y r bm *» y * 1 ■ y * * < fc.o fci55«<D«*i=B 

[0 0 4 8] ±EB. KT?ff5C i 0>E$r*&3I* 

^y-UBftfUt^ii:, XU*K**#THRE 
AD<fcLT. L* 2 x#TH READ<llC#S]L. CCD 
T=.fW&*>mC< 2x#THREAD«|:»at4« 
^Lt, ±iT^&#THREADi0'<7*ft«, - 

[0049] m 161*. ld lT ^«(DHJSMT;u=ry 
XA^t^-eft^o ldlT ##<dt;u zjyxAi*. 

It^iU— ^>LDLtLTS8a**L*. 

s&*>£»(*. «raa)Lu»»a)a^i:i^«|-efe4o 

[0 0 5 0] nwid*<8<0»*l::i** (2 0) 0)ffl5 
ESldJ #S*155o -t L"C. LT i (ist 
+ 8 : k. ist : ist +7) £JE$rf*o -<&<t£. L 
T i (ist : ist +7. ist : ist +7) 0>± = ftgfl# 
ICDLT 36<AoTl^O>"C. (DLT ) -1£&&\C>AM+ 

[005 1] (2 0) 4>M<DI FXt, nwidA<8-eft 
l*4r«K**ifc*6l-te* nwid/2**rfcttnwid4: It 
f>LD L*a?t;ttiL. *fff*. L 
T i (ist : ist +nwid/2-1. ist +nwid/2 : is 
t +nwid-1) IcdlT *:□£—-*-*. Dl*. LT i 
(ist : ist +nwid/2-1. ist : ist +nwid/2- 
1) Ottftmmv&V. LI*. LT i (ist +nwid/ 
2 : ist +nwid- 1 . ist : ist +nwid/2- 1 ) ^fc 

y. zajLseitSo 

[0052] fit, LT i (ist +nwid/2 : k. is 
t +nwid/2 : ist +nwid- 1 ) £EJT*"£o "Ta*> 

LT i (ist +nwid/2 : k. ist +nwid/2 : ist + 
nwid-1) =LT i (ist +nwid/2 : k. ist +nwid 
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/2 : ist +nwid- 1) -LTi (ist +nwid/2 : 
k. ist : ist +nwid- 1 ) x LT i (ist : ist +nw 
id/2 - 1 . i st +nwid/2 : ist +nwid- 1 ) 

[0 0 5 3];^:, ldJ ^SOT^u-^vldl 

sic* **>qEynx*^iMmjMM> 

■JXAI*. CD-ROM^WiftEflkaf*^ /\— Kf 

[0 0 5 4] za>*aft»*. **K<7>ll!6fl2®<07;u 
(##£») 

■ LU5M*a>£»*. 1) . 2) % »3UX4-^S 

1) P. AMESTOY, M.DAYDE, and I. DUFF, "Use of c 
omputat i ona I kerne I s i n the so I ut i on of f u 1 1 
and sparse I inear equations" , M. COSNARD, Y. R 
OBERT, Q.OUINTON, and M.RAYNAL, PARALLELS; D I STR 
IBUTED ALGORITHMS, North-Holland, 1989, pp. 13- 
19 

2 ) G. H. Go I ub, C. F. van Loan, "Matr i x Computat 
ions" , second edition, The Johns Hopkins Un 
iversity Press, 1989 

• B*H(D*:«Tftt. jaTOtOlZLUfltfXtfLDL 

■ rjaflipwj « jes*. &&&&&& 

P^^Ort. »ft^P^*£»fl:7P^£T?ftlv]\fTM 
y*ij««ic»tt-r*»ttx^^?i:. ?P-fcr 
nTZ>zt\z±v. «ao^p-fe^it^^3.-;u-e. » 
t" ^-et# b <xfc /hff n ^ p ? £ <a?li$r f§s £ m o r % 
tt-4it«fr«»a*a*i«*B««i3*a*i*4^py 



[0 0 5 5] (f*E2) ttfrM*)tf±» ff«0>LU5MB 
T*fc4 - £ £4#® fc-T 4ttE 1 [cSH«0)I5«i«<* o 

(^IB3) 8frE«!S<D^P^^-9" : tvj.--;u^-r4/h 
fi^J^P^^OT^r-^^b^ix^^tfTtt^ KD0Mi£ 

O^Piz^^^^-^^^jS^^^U^Iit^fcLN-C. 
t f 4ttE 2 lcEtt<Z>ES«Et*o 

[0 0 5 6] (#E4) fflELU#Sli. HrEflM©* 

e»n»M*r;u=f u xakj: o -ciBsMrJiJoEtr* 

^^P-b^-tr^vzL-^U-C-LU^S-r^CirlC^y. tt 
ftE2lzEK0)Efa«ffo 

[0 0 5 7] (#E5) KfrMlftf** 

WE6) «9En l/X+-Mft^l>l«r^ L/X*- 

* its t -r 4f*E 2 izEis<7)Eiij«i*o 

[0 0 5 8] (f*E7) BrESSr^^^^lzfcixr. S 
•r-r^S^frW ^P ^ * <&HS*T?]ffl5# * iESa© 7 
P-b ^-tr^E Vi— Jl(D&<D 2fe<D®L<Dyn V*? 
U K^f»J^*tfc = ft^JaJ^O^P^^* 2 0-foijfi 
H^io-tex^ ^^p-tz^-lJ-^EvzL— ;u(DP— ^EU 

t**«i:-r*#E5lzE«a)Ett«#. 
[0 0 5 9] (#E8) ISjaO^P-iz^-tr^Evi-^S 

O X I vhfr^J ^ p ^ ^ £ * KSSfeO ^nt^^yi 

va- ;utj, »^^p^^^7i:S(3;1IS:-r4;1IS:x-r^ 
ft. mm&XTV?xm*tit^'n&\7nv2<»8iU 

[0 0 6 0] (f*E9) S^CD^P-tr^-tr^vzL-^^ 

^to^> ; Eyax*^^Jtt®:*ilcfcL^T. ffjit* 

^T^lj ^ P ^7 fC» It 4 37 P ^ 1 1b^K it . K/J^tT^J ? P 



(8) 



$$[182002-1 63 246 



[006 1] 

u X -r a) & £ ftw <*> a m ysA & *t a . 

[01] *#*^US^*9tt*JtW«a>>\-K*xT 
[0 2] **W0)*16»Blzfl£ofcLU»»a>*WAa 

[0 3] **B^(DHJfi^BI-ttofcLU^ScDM^JAS 
O)«E^:*SiBJ-r^)0 (*<7)2) T?*«o 

[0 4] *^J6^®(DLU^»^AS<D35n*^-r«tBS 
?P— ^-V — hTrfc&<> 

[0 5] MlkMa)LU«flia)«tt*J:UIMicHM 
-r^>0 (-ta>i) 

[06] 2Mm«a0Lu^a>»ft**yiMm=m 
-r&0 (^2) r*fe*o 

[0 7] *36K^»o)Lu»»a)***cfcypfflicKw 



[0 8] 4:iiik»fioLu««o>«tt«<feumai=Km 

[0 9] **j6^ffia>Lu»»a>^**y»fflicjttP 
-r*B c^cds) t?ft^ 0 

[01 o] *se«E^ffla)Lu»»(z>^a*<fcyBfflicift 

^-r^>0 (*a>6) -efe^o 
[011] jE{I#ffcfr?iJ<Di«£lc=i U>X*-#M£*t5 

A«a>*«:*lMfl"r*Bi (-tan) -efe*o 

[01 2] ]Eflia»ffW<01i$lcaUX+-»«Sff5 
[01 3] IEtettft^J<Di#£iZ=l UX^-^^-^a 

Aaa>*)fc*Kwr4H (*a>3) -cfc^o 

[01 4] t«3 l/Xt-iMBOTi^UXA* A y B 

[015] x»3 i/x^-^tj^u xas * tie 
iz&BXi-zm (^2) -cfc*. 
i e] K»aux*-#»a>TJU^"JXA**y» 

10-1-10-n ^P-fe^it 

1 1-1-1 1 - n ^*'Jt5a- ;U 

1 2 ffiSS^S 

1 3-1—1 3-n 2**-V75/a/*U 



[01] 



[02] 



IH H-2 



11-n 



12 



10-2 



in 









Mi****) 



10-n 









\>uu 
ll\ 


: u ; 
' — j — " — ! — - 

Ul | U2 | U3 




LI 


Ct 




L2 


C2 




L3 


C3 1 
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[El3] 



4] 



D 




D 




D 


L1 




LZ 




L3 


Tl 


T2 


T3 




[05] 



•4 



448 



448 





U 




Dt 


UI | 02 J U3 


U4 


LI 


ct 




L2 


C2 


L3 


C3 


L4 


C4 



f_ START 



I 



S1 



32 




I 



86 



C 



I 



END 



3 



[06] 



K*1 *U>*ft2 Xl^t3 3U»pi« 



01 




Dl 




Dl 




D1 


L1 




L2 




LS 




L4 
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[07] 



[08] 




*3Ot»»<0LjU#IB(D;5i££ 
iyB«=a*-*-*BB(*©3T 



OO 5-1. tblc, 

TMP»0.0 00:|j«so 
DO j-l. hnc 

LTQ. J)}. QT . tmp;then 

TMP-ABSCLTO. 0) 

JJ-i 

END1F 
ENDDO 



<i> 




ifOi. arr, o then 

DO k-1. ibtks 
TMPX«LT(1. 10 
LTCl.k)«LT<B.k) 
LT(fl,k)»TlWPX 
ENDDO 
END IP 

DO K-l+1. IbBcs 

LT(l.k>=LT<l.k)/LT{l. i> 
ENDDO 

DO k~l+1. ibBcs 
DO 1«l + T,i»n« 

LTO. k)»LT(l. W-LTO. •} XLTO, k> 
ENDDO 
ENDDO 

ENDDO 



(2) 



<3> 



DO 1=V 
DO 

WLAM LTia I)). ST . TM Pi THEN 
TMP«AM(I-Tia, B) 



p*«p«C(*T»«?*A0)~| 
((PTMRHAJJtt- Xl/->*>»4. 4XI/3+* 



] 
] 
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(B) 



IF(WTHRKM>. «Q. 1 > 

J ebO; GPM9T «0 
DO k(*1.4 



F{*-mRCAD, BO, «JPWOT>TMEM 

iFTj.arr, vtwn 

DO bta 

ROWOO-LTItt. hO 
ENOOO 
END T 



(0> 



C7) 



F(GWOT. EQ, OTTMBN 

rxii. arr. stub* 

DO l-t. Mft 

t up*- una ik) 

LTiaW«lTi<J. *) 
UTlOj.faO- 



. f juh»JU*rt©i«f*9fc 1 



END IF 



CTPTHREMX ED, G PIVOT) TMt* 
DO M* 

una *o«una w 
unci. m-rowiw 



ELBE 

DO k>V Un 



DO k-i+T.tek*. 

LTIU W-LTiU tO/LTtl. D 
ENDDO 

OO k-i+I.B*- 
DO J-H-1, tan* 

una w-L-nn. k>-una o ««-Tm w 

ENDDO 
ENDDO 



] 
] 
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£101 



CNUOO 



[09] 

tfcuH T*12 (^osi 







384 

v, 

r 

384 

\ 


01 


U1 | U2 | U3 U4 


Lt 


C \ 


L2 


C2 


L3 


C3 


L4 


C4 



m i o] 



•ubroutln* UJCLTl, k. Uks, W. mrfd) 
(LTI»#X^PTD1 +U*Mft.- 

IFCnwfd, «q. a}TWt(jaafiS»^ 
LTlCat:k. i»t i*l+nwkl-1)*!k#M=lJU*MfT«. 

r zc-c, (4>-(io)«r5. I 



cai UUCLTl, k. ft**, tat nrwidX2) 
call TRBC ) 

LT1 (bW : tat+rn«*d/2- 1 , |*H-rwW/2:l«t+nwk0 
tllti. LTUI*:l»H-fmkl/2-1, Ufc «t+mrW/2- 1 ) 

c^t MM( ) 

CLTia«H-nwkl^£:lw l«t-Hnwld/e : l^+m»ld> 
■L71< fc d; I nwi d /2:K. l*H-nwfcl/2:t«t+rrwkl> 
-LT1(b*+nwkl/2:k. l*:!it+Twrid/2-1)>« 
LTi(i*:lrt+nwidX2-1. tafc+nwW/2:tal+nwld> 

call LUCLT1. k lb**. bH-ftwd/2, nwid/2 
•nd if 
ratum 

•nd .ubnoudn* 



] 
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1 1] 



[1121 




[113] 



■ 









SI 



T /* 









I I 










III' 1 




1 




o 











l-M 1 



a » 











uu 




-H \ 
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4#B32 0 0 2- 1 6 3 2 4 



1 5] 



DL2I 




[Hi 6] 



cufaroutfea LTCXLTI. k. IWWi. lat n«4d) 
BECnwid. EQ, 83 THEN <«lJ8tf»'M 
DOi-iBt lat+7 
D04»l+1,l«t+7 

LTKI,i>=LH(i Q 
LTIO. 0-LTia 0/»-T4{l, !) 
ENDOO 

DO Jy—l+1. 1«+7 
DO iz^lx, k«t+7 

LT1 Q*. Jy) ™LTi Ox. jrf -L.T1 (j* •> W LTiti. far) 



C20J 



EMDDO 



J 



CLTI(U7iCirt+B;k, W:!»t+7)tWlf ^ 

call LDL<LT1. k. Iblkm, lat. nwW/2) 

.i n/i«t l«M.n-trf^g— I l«t+m«>dX2 :l«t+rrwid- 1 ) 
!=DL T eat;-«- (DttLTi(lrt:tet+rmid/2-1. tat:tat+nwid/2-1) 



LTKIa*+nwM/2:k, l«t+nwU/2:i«H-™rid-1) 

=LTI CM: k. n*kl/e:iat+ wM- 1 ) - 

L71(i»l+nii*iy2:k. Ia»:i«t+rwiikl-1) « 

CAUL LDL (LT1. k. ibika, iat+nwrid/2. 



] 



RETURN 
END 
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